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CN . Abstract 
> . 

Q , Charmonium sum rules are analyzed with the primary goal to obtain the restrictions 

' on the value of dimension 4 gluon condensate. The moments M„((5^) of the polarization 

operator of vector charm currents are calculated and compared with experimental data. 
! The 3-loop (a^) perturbative corrections, the gluon condensate contribution with ag 

corrections and dimension 6 operator contribution are accounted. It is shown that 
the sum rules for the moments do not work at = 0, where the perturbation series 
. diverges and contribution is large. The domain in the plane {n,Q'^), where the 

I sum rules are legitimate, is found. Strong correlation of the values of gluon condensate 

' and MS charm quark mass quark is determined. The absolute limits are found to 

. be for gluon condensate (^G'^) = 0.009 it 0.007 GeV^ and for charm quark mass 

g ; rh{fh) = 1.275 ± 0.015 GeV hi MS scheme. 



^ ; 1 Introduction 



PACS: 14.65.Dw, 11.55.Hx, 12.38.Bx 



It is well known, that QCD vacuum generates various quark and gluon condensates, the 
vacuum expectation values of quark and gluon fields of nonperturbative origin. Among 
them the gluon condensate (^C^j.G^^), where is the gluon field strength tensor and 

2 

as = ^ is the running QCD coupling constant, plays a special role. The existence of the 
gluon condensate in QCD was first demonstrated by Shifman, Vainstein and Zakharov 
Its special role is caused by few reasons. First, it has the lowest dimension d = A among 
gluon condensates, as well as any other condensates conserving chirality. For this reason 
the gluon condensate is the most important one in determination of hadronic properties by 
QCD sum rules, if chirality conserving amplitudes are considered (e.g. in case of meson mass 
determination). Second, the value of the gluon condensate is directly related to the vacuum 
energy density e. As was shown in |]l[ 

^ = (1) 

where is Gell-Mann-Low /3-function. Therefore, the sign and magnitude of (^G^) are 

very important for theoretical description of QCD vacuum and for construction of hadron 
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models (e.g. the bag model). Third, in some models the numerical value of the gluon 
condensate is usually used as normalization scale, which fixes the model parameters. For 
example, in the instanton model it is required, that this value is reproduced by the model. 
The numerical value of the gluon condensate 

f = 0.012 GeV^ (2) 

has been found in |I| from charmonium sum rules. (This value is often referred to as the 
standard or SVZ value.) Later there were many attempts to determine the gluon condensate 
by considering various processes within various approaches. In some of them the value (0) 
(or ones, by a factor of 1.5 higher) was confirmed in others it was claimed, that the 



actual value of the gluon condensate is by a factor 2-5 higher than (H) [[^- [[14 |. 

From today's point of view the calculations performed in |^ have a serious drawback. 
Only the first order (NLO) perturbative correction was accounted in 0] and it was taken 
rather low value of a^, later not confirmed by the experimental data. (It was assumed, that 
QCD parameter A'-^^ ^ 100 MeV and ^^(mc) ~ 0.2, today's values are essentially higher.) 
The contribution of the next, dimension 6, operator was neglected, so the convergence 
of the operator product expansion was not tested. In charmonium sum rules the moments 
MniQ"^) of the polarization function n(g^), = —Q^ were calculated at the point = 0. 
It was shown in [jl4|, that the higher order terms of the operator product expansion (OPE), 
namely the contributions of and G^ operators are of importance at = 0. The results 
of calculations of the second order (NNLO) perturbative corrections to n(g^) as well as a^- 
correction to the gluon condensate are available now. They demonstrate, that both of them 
as a rule are large and by no means can be neglected in the sum rules for the moments at 

= 0. Finally, the experimental data shifted significantly in comparison with ones, used 
in 0. 

Later the charmonium sum rules were considered at the NLO level in for > 
and their analysis basically confirmed the results of There are recent publications [^, 
T6[] , p!7[] where the charmonium as well as bottomonium sum rules were analyzed at = 



with perturbative corrections in order to extract the charm and bottom quark masses 
in various schemes. The condensate is usually taken to be or some another fixed value. 
However, the charm mass and the condensate values are entangled in the sum rules. This 
can be easily understood for large Q^, where the mass and condensate corrections to the 
polarization operator behave as some series in negative powers of Q^, and one may eliminate 
the condensate contribution to a great extent by slightly changing the quark mass. Vice 
versa, different condensate values may vary the charm quark mass within few percents. 

The condensate could be also determined from other sum rules, which do not involve 
the charm quark mass, but the accuracy usually appears to be rather low for this purpose. 
In particular, precise analysis of e~^e~ data [0 lead only to rather weak restrictions on the 
gluon condensate. In ref |T^ the thorough analysis of hadronic r-decay structure functions 
was performed and the restriction (^G^) = 0.006 ± 0.012 GeV^ was found. This value, 
however, does not exclude zero value of the condensate. 

For all these reasons a reconsideration of the problem is necessary. The charmonium sum 
rules on the next level of precision in comparison with |^ is presented below. In Section 2 
general outline of the method is given and the experimental input data for the sum rules are 
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presented. In Section 3 the method of calculation of the perturbative part of the moments 
is exposed with the references to the sources, we used in the calculations. Section 4 presents 
the gluon condensate contribution with a^-corrections in the form, convenient for numerical 
evaluation of the moments for nonzero . In section 5 perturbative and operator product 
expansion of the moments is considered. It is argued, that the choice of pole charm quark 
mass as a mass parameter is not suitable, since in this case the higher order in terms 
overwhelm the lower ones and the a<j-series are divergent. It is proposed to get rid of this 
problem by using the MS mass as the mass parameter. In what follows the MS charm quark 
mass m[m) at the renormalization point equal to the mass itself is used. The formulae for 
moments Mn{Q'^\ expressed through the MS mass, are given and the domain in in^Q'^^ 
plane was found by direct calculation, where the perturbative series are well convergent. 
In Section 6 the calculation of m(m) and gluon condensate is presented. In Section 7 the 
sensitivity of the results to the operator contribution is tested. Section 8 is devoted to 
the discussion of the attempts to sum up the Coulomb-like corrections. Section 9 contains 
the conclusion. 



2 Experimental current correlator 

Consider the 2-point correlator of the vector charm currents: 

dxe"^"" {TJ^{x)J,{Q) ) = ( q^q, - g^,q^ ) Tl{q^) , = c^^c (3) 



The polarization function n(g^) can be reconstructed by its imaginary part with the help of 
the dispersion relation: 

R^{s) = 47r Im n(. + .0) , U{q') = T • (4) 

4vr^ i4m2 s{s - q^> 

We shall use notation Rc not to confuse with frequently used notation R for the imaginary 
part of the electromagnetic current correlator, R{s) = 3(5ji?/(s), the normalization 
Rc{oo) = 1 in the parton model. In the narrow-width approximation Rc{s) can be repre- 
sented as the sum of the resonance (5-functions: 

3 TT 

^^(^) = n2 2 ( \ y'm^r^^ee5(s-mj) (5) 

<5cttcm(s) ^ 

where Qc = 2/3 is electric charge of c-quark, a^rnis) is the running electromagnetic coupling: 
, a(0) 



Aa{s) ' 



Aa{s) = -47ra(0)ne„,(-s) = Aaiep(s) + Aahad(s) . (6) 



Here a(0) = 1/137.04 is the fine structure constant, nem(s) is the correlator of electromag- 
netic currents J^™ = J^iQii'ilfi^'i defined in the same way as (j^). As usual, the leptonic 
contribution to nem(s) is found by the perturbation theory, while the hadronic contribution 
has to be determined by numerical integration of experimental e'^e~ (or r-decay) data. Since 
Ciem{s) wcakly changes from one resonance to another, we fix it at s = 'Tijy^ from now on: 

Aaiep(mj/^) = 0.016 , Aahad("^j/^) = 0.009 , aem{m]/^) = 1/133.6 
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Figure 1: Rc{s) in the region of high resonances, determined from BES data pT 



There are 6 vector charmonium states with J^*" = 1 [2^: 



notation 


mass, MeV 


full width, MeV 


r^^ee, keV 


jms) 


3096.87 ±0.04 


0.087 ±0.005 


5.26 ±0.37 


V(25) 


3685.96 ±0.09 


0.300 ±0.025 


2.19 ±0.15 


V'(3770) 


3769.9 ±2.5 


23.6 ±2.7 


0.26 ±0.04 


^(4040) 


4040 ± 10 


52± 10 


0.75 ±0.15 


^(4160) 


4159 ± 20 


78 ±20 


0.77 ±0.23 


V'(4415) 


4415 ± 6 


43± 15 


0.47 ±0.10 



The first two resonances, J/tp and tp{2S), are sufficiently narrow and their contribution to 
Rc{s) can be well parametrized by the 5-functions (|). 

But the next resonances, especially the last three ones, are rather wide and the narrow- 
width approximation for them could be inaccurate. Here it is better to use Rc{s), extracted 
from e~^e~ — > hadrons branching ratio R{s) = 3Q'jRf{s), which is measured experimen- 
tally in wide range of s. Precise data on R{s) in the region of high charmonium states were 
obtained recently by BES collaboration [^]. In order to extract Rc{s) from these data, one 
has to subtract the contribution of the light quarks from R{s). We suppose, that it is well 
described by the perturbative QCD, which gives 2.16. The result for Rc{s) is shown in Fig 
|l]. Above the last resonance Rc{s) is getting close to 1, the parton model prediction. 

Now we summarize the following experimental input for Rc{s), which will be used in our 
calculations: 

s < Si = (3.7 GeV)^ : (5-functions from J/ip and ip{2S) according to eq (5) 
s^<s<S2 = (4.6 GeV)2 : BES data, see Fig |I] (7) 
S2 < s : continuum, Rc(s) = 1 



One could include the a^-correction in the continuum region, but they will not be essential 
in what follows. 



4 



In order to suppress the contribution of the high energy states, one considers the deriva- 
tives of the polarization function in euchdean region = —Q"^ < 0, the so-called moments: 

The experimental values are calculated according to (0): 



[s + Q 



2V«+1 



M (O^) = V ^^^^^^^ I r + 1 (q] 

4 aL ^ K + Q')"^' ^ + ^')"^' ^(^2 + ^ ' 

The squared error of the moments (|^) is computed as the sum of the squared errors of each 
term. 

The lowest state J/ ip gives maximal contribution to the moments due to the largest width 
Tj/^^ee, which itself has the error 7%. This error can be eliminated to a great extent, if one 
considers the ratio of two moments, which in general case can be written in the following 
form: 

^<"""^^'?^>^tS!^<"5/* + o^r-"'iif . (10) 

where D12 denote the higher state contribution to the moments divided by the J/ip 
contribution. Then the error of this ratio is calculated by usual rules: 

where the mass errors are neglected. If D12 ^ 1, the relative error of the ratio is much 
smaller than the relative errors of the moments itself. This fact has been utilized in many 
papers on charmonium sum rules and will be used here. 

In our calculations we shall always use sufficiently high n (n > 8), so that the last term 
in (^, which comes from continuum, is small compared to the resonance contribution and 
the uncertainty introduced by this term is negligible. Moreover, the difference between the 
narrow width approximation for the high resonances (above il){2S)) as given by (^, and 
their representation by (H) is small and well below the quoted errors. 



3 Theoretical R{s) 

At first one defines the running QCD coupling a(/i^) = as(/i^)/7r as a solution of the renor- 
malization group equation: 

Then the functions R^"'\s , fi^) are defined as the coefficients in the a^-expansion: 

R,{s) = 5^i?(")(s,/i2)a"(/i2) (13) 

n>0 
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Since Rc{s) is the physical quantity, it does not depend on the scale /i^, although each term 
in ( ]T3D may be /i^ dependent. 

It is easier to represent the results in terms of the pole quark mass m and the velocity 
V = -^/l — 4m^/ s. The first two terms in the expansion (|l^) do not depend on /j^. The 
leading term R^'^^ was calculated in the next to leading R^-^^ in [p^ : 



_(5_3^;2) + 2^(3-^2)(ln ^ 



l + v 



2Li2 



1 - t; 

l + v 



Lin 



(14) 



where Li2(x) = Yl'^=i ^ the dilogarithm function. The function R^"^^ is usually decomposed 
into four gauge invariant terms: 



CpR^2^ + CaCpR^^\ + CpTriiR^"^^ + CpTR^P 



(15) 



where Ca = 3, = 4/3, T = 1/2 are group factors, rii = Uf — 1 is the number of light 

(2) 

quarks. The function R^ comes from the diagram with two quark loops: one loop with 
massive quark, which couples to the vector currents, and another massless quark loop (the 
so-called double bubble diagram). It was originally found in and in our normalization 
takes the form: 



(2) 



4 4s 



5 

12 



(16) 



where the function 5^'^^ is given by equation (B.3) in ref [^]. The function i?^^ comes from 



the similar double bubble diagram with equal quark masses and has the form |26|: 



R 



(2) 



In 4 

4 m'' 



(17) 



where is given by equation (12) in ref [^]. The function comes from the 4-particle 
cut and vanishes for s < 16 m^. It is represented as the double integral (13) in ref [^] which 
can be computed numerically. However for s > 16 the total function Rp^ can be well 
approximated by its high energy asymptotic: 



R 



(2) 



In^ 

4 s 



+ Cs 



11 



13 m 



In numerical calculations we take all the terms up to m^'^/s^, extracted from [^]. The 



(2) (2) 

functions R\ and -R)vyi generated by the diagrams with single quark loop and various 

(2) (2) 

gluon exchanges, R\ is abelian part while -R)v^ contains purely nonabelian contributions. 
They are not known analytically. We will use the approximations, given by equations (65), 
(66) in ref |25| (divided by 3 in our conventions) which reproduce all known asymptotics and 
Pade approximations with high accuracy. 
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4 Condensate contribution 

The contribution of the dimension 4 gluon condensate (aG^) = (^CJ^^G^^,) to the polariza- 
tion function of massive quarks has the form: 



{aG^) 



(4m2)2 

The leading order function was found in |]1| 

1 



(1)/ 



where 



Am? 



12 z^v^ 



3 2z - 1 , V 
m 



zv 



v + l 



z + 



(19) 



where v = \/\ — 1/ z. For this function the following dispersion-like relation can be written: 



1 



dz' 



12 ii z'-'v' 



4 iz' - z) 



+ 



[z' - 2; 



(20) 



This representation is convenient for an evaluation of various transformations of the polar- 
ization function n(s), in particular, the moments. 

The next-to-leading order function f^^^ was explicitly found in [0. One could differen- 
tiate it n times to obtain the moments for arbitrary Q^. However, we prefer to construct 
the dispersion integral similar to (^). The function f^'^^z) has a cut from z = 1 to oo and 
behaves as at z — > 1. Integrating f^^\z')/[z' — z) by z' along the contour around the 
cut, one obtains the following representation: 



1 

+ 



lm/W(z' + iO) 



8633 



z' — z 



i=i ^ 



65 



1152 I- z 



17 

+ — In (8e) 
6912 36 ^ ' 



.-1/2 



65 



'1/2 



384 (1 - z)'^ 



where e — and 



/i 



17 
384 



/2 



413 
6912 



/3 



197 
2304 



(21) 



(22) 



The imaginary part is: 



Im/W(;z + iO) 



TT 



962;^f ^ 



+ Pl{z){l-z) [ 21nw + ^ln(4z) 



zv 



l + v 



(23) 

where the polynomials P^{z) are given in the Table 1 of ref |jl2|. It behaves as f ^ at 2; ^ 1, 
so the integral in ([21| ) is divergent in the limit e — > 0. We decompose it into 3 parts 



hmf^'\z + ^0) = Fr{z) + F^{z) + ^ F^'(^) 



TT 



(24) 
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in such way, that each function Fi{z) behaves as w ^ at 2; — 1 and has appropriate asymptotic 
at 2; — > 00. In particular we choose 



FAz) 



1 



96 z^v 



1627 893 



36 



+ {--22z 



9 
140 



98 , o 
-z - —z^ + 368;z^ 



3208 



'-z^ + —^z^ - 3232^^ + 2208z^ ) — In ■ 

3 J zv 1 + V 



1 . l-v 



+ ( 682 + ^z^ - ^2=^ + I4722M ( 21nt; + - In (4^; 
9 3 / \ 2 



96 z^i; 
136 



23 559 



-z + 



272 



197 947 \ 1 1-t; 
1 z — In 

8 12 zv l + v 



FAz) 



+ ^z"^ l2\nv + ^ln(4z) 



67 197 , l-v 
in 



96 z^t 



6 12 zv l + v 



(25) 



Then one may integrate (^) by parts twice and all singular in e term cancel. Eventually we 
obtain following representation for the function f^^^: 



i=l 



[{i-zy 



+ 



[z' - z) 



dz' 



(26) 



It will be used to compute the moments numerically. 



5 Moments in MS scheme 

For definiteness let us choose the scale /z^ = rn? in 
the moments (|): 



and write down the Os-expansion of 

(27) 



A;>0 



k>0 

The perturbative coefficient functions are 



B!h\s, m^) ds 

The leading order can be expressed via Gauss hypergeometric function: 



ii^°)(g^ 



1 3v/? (n + l)r(n) /n,n + 2 
(4^2)n 4 r(n + 5/2) V^ + 5/2 



4^2 



(2^ 



(29) 



The higher order functions M^^^ and M^"^^ are computed numerically by (28). (Notice, that 
the analytical expression for Mn\o) has been found in and the ffist 7 moments Mn\o) 
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can be determined from the low energy expansion of the polarization function 11 (s) available 
in 



The leading order contribution of the gluon condensate is easily obtained from (pO]): 



TT 



{n + l)r(n + 4) / n + 2, n + A 
(4m2)«+2 "6 r(r2 + 7/2) ^ \ n + 7/2 



4m2 



(30) 



The next-to-leading condensate correction can be computed numerically with the help of the 
integral representation, obtained from (p§): 



Mf'i)(g^) 



Aw 



{Am? 



n+2 



i=l 



T{n + i) 
T{n + l)TU) 



Fi(z) 



(z + y) 



dz 



(31) 



where y = Q'^/{Ani?)^ the constants fi and the functions Fi{z) are given in (|2^) and (p5|). 

The pole quark mass m is the most natural choice, since it is the physical invariant. 
However in the pole scheme the perturbative corrections to the moments are huge. For 
instance, at the typical point, which will be used later in our analysis, one gets: 



n 



10 , 



Am 



13.836 



M(o) 



193.33 



13.791 (32) 



Since in the domain of interest a ~ 0.1, this is an indication, that the series ( pTj) is divergent. 
The situation is even worse for = (see |]l5l). It is almost impossible to choose an 
informative region in the {n, Q"^) plane where the perturbative corrections in the pole mass 
scheme are tolerable and the continuum as well as {G'^) contributions are suppressed enough 
on the other hand. 

The traditional solution to this problem is the mass redefinition. In particular, in the 
most popular MS scheme the mass corrections are known to be significantly smaller. In MS 
conventions the mass fh depends on the scale iJ? according to the RG equation: 



m{ii^) = m{fil) exp ( / '^wr^'^a 



7m (a) = y^7nQ 

n>0 



n+1 



(33) 



where 7^ is the mass anomalous dimension. In what follows we shall choose the most natural 
mass scale fi^ = m? and will denote fh{m?) as simply m. 

There is a perturbative relation between the pole mass m and MS one fh: 



m 



m^ 



1 + ^ ir„a"(m2) 



(34) 



n>l 



The 2-loop factor was found, in particular, in while the 3-loop factor was recently 
calculated in |p9| : 



K2 



8 
3 

28.6646- 2.0828 = 22.4162 

417.039 - 56.0871 n, + l. 3054 = 260.526 



(35) 
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-1 



5 10 15 20 5 10 15 

Figure 2: Ratio mI^^ / m'^^ (left) and mI^^ / M^^^ (right) in the plane (n, Q^). 



20 



We put = 3 in the last column. The series (34) also looks divergent at the charm scale. 
(Notice, that the authors of |I| used another mass convention, although numerically close to 
MS scheme at the NLO level: the coefficient Ki was equal to 4 In 2 there.) 

Nevertheless let us assume for a moment, that is small, take advantage of (P^,pH) and 
express the moments (P) in terms of the mass m: 



k>0 



(36) 



k>0 



As follows from the definition 



and dimensional consideration 



M^i) - 7^1 (n - d/2) M^°) + K^{n + 1) Q^M^^Ji 



M^2) ^Ki{n- d/2) M^^) + Ki{n + 1) Q'^mI^I, 
2 



+ (n - d/2) 



-^{n + l-d/2)-K2 



Mf'°)(g^) 

Mf'i)(Q^) 



+ (n + 1) [7^2 - Kl{n + 1 - d/2)] Q^m'^I + ^^^(n + l)(n + 2) Q^M^^ 



n 

Mf -Ki{n + 2- d/2) + iri(n + 1) Q^M, 



(G,0) 
n+l 



(37) 



where d is the dimension of the polarization function n((5^) (c? = for vector currents), all 
Mn ■* in the rhs are computed with MS mass fh. 

The moment corrections M'^^'> are much smaller than M*^^) in the pole scheme. In partic- 
ular, at the same point, which was considered in (|3^), we have now: 



n 



10 , = 4m 



M(o) 



0.045 



1.136 



-1.673 (38) 
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5 10 15 20 

Figure 3: Ratio Mf'^^/Mi^'^'^ in the plane (n, Q^) 



This smallness of corrections as compared to the pole scheme is observed for almost all n and 
Q2. The ratios M^^V^n°^ and Mr / M^^^ are shown in Fig | and the ratio mI^'^V^*^^'"^ 
in the Fig ^ for n = ... 20 and Q^/(4m^) = —1 ... 3. The perturbative expansion in MS- 
scheme obviously does not work in the area of high n and low Q^, marked with dark. (The 
detailed data are presented in the Tables 1,2,3 in the Appendix.) 



Now we can argue, why the expression (|37|) for the moments is legitimate, despite that 
the series (|3^), relating the pole mass m and MS mass m, is divergent at the coupling as 
taken on the charm mass scale. If as is small enough, eq (|37|) is correct. In this case the same 
values of can be obtained by the procedure, when MS mass renormalization is performed 
directly in the diagrams, without all the concept of the pole mass. If the pole mass concept 
is not used, the relations (|3^ , |35D are irrelevant. These relations demonstrate only, that the 
pole mass is an ill defined object in case of charm. The check of selfconsistency of MA 
moments is the convergence of the series (pGl). 

If one takes the QCD coupling at some another scale q;s(/U^), the function M*^^^ must be 
replaced by: 



(39) 



so that the series ( ^6|) is /i^-independent at the order a]. 



6 Determination of charm quark mass and gluon con- 
densate from data 

Theoretical moments depend on 3 parameters: charm quark mass, QCD coupling constant 
and gluon condensate. The QCD coupling is universal quantity and can be taken from 
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■( 22 ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' I ... I ... I ... I ... I ... I . 

6 8 10 12 14 16 18 20 -0.03 -0.02 -0.01 0.01 0.02 0.03 



Figure 4: a): MS mass found from experimental moments M„(Q^) for different n and Q\ 
determined by the equation Mn\Q'^ = for different values of the gluon condensate. The 
shaded area shows the experimental error for <^^G^) = 0, for nonzero condensates only the 
central lines are shown, b): fn{fh?) in GeV vs l^^G'^'^ in GeV^ determined from n = 10 and 
Q2 = 0.98 X Am?. The is taken at the scale (pj). 



other experiments. In particular, as boundary condition in the RG equation ([1^) we put: 

as{ml) = 0.330 ± 0.025 , = 1.777 GeV (40) 



found from hadronic r-decay analysis [|T^ at the r-mass in agreement with other data ||20[ . 



Another question is the choice of the scale /i^, at which as should be taken. Since the 
higher order perturbative corrections are not known, the moments M„((5^) will depend on 
this scale. In the massless limit the most natural choice is fi^ = Q^. On the other hand 
for massive quarks and = the scale is usually taken /i^ ~ m^. So we choose the 
interpolation formula: 

= + rn" (41) 

At this scale as is smaller than at /i^ = for the price of larger Mn^ according to (|39|). 
(Notice, that in the Tables in the Appendix as well as in the Fig |^ the ratio M^'^^ M^^^ is 
given at the scale /i^ = m^.) Sometimes we will vary the coefficient before (O) to test 
the dependence of the results on the scale. 

The sum rules for low order moments M„((5^), n < 3 cannot be used because of large 
contribution of high excited states and continuum as well as large a^ corrections (see the 
Tables in Appendix), especially at = 0. As the Fig | demonstrates, at n > 4 the 
correction to the gluon condensate is large at = 0. The (G^) condensate contribution is 
also large (see below), which demonstrates, that the operator product expansion is divergent 
here. For these reasons we will avoid using the sum rules at small Q^. 

As the Fig |^ shows, the first correction to the moments Mn\Q'^) vanishes along the 
diagonal line, approximately parametrized by the equation Q'^/ (4m^) = n/5 — 1. The second- 
order correction M^"^^ and the correction to the condensate contribution M^*^'^) are also 
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-0.01 -0.005 0.005 0.01 0.015 -0.005 0.005 0.01 0.015 0.02 



Figure 5: MS mass versus gluon condensate obtained from different points on (n, Q^) plane. 
"Horizontal" bands obtained from the moments (10, 4m^) and (15, 8m^), "vertical" bands 
obtained from the ratio of the moments M10/M12 (left), M15/M17 (right) for few different 
choices of as(/i^). 



small along this diagonal for n > 5. Now lest us compare the theoretical moments with 
experimental value (P) at different points on this diagonal. If the condensate is fixed, then 
one can numerically solve this equation in order to find the MS mass. The result is shown in 
Fig The values n < 5 are not reliable, since the Os-correction to the condensate exceeds 
-50% here. 

The lines in Fig^a are almost horizontal, if the condensate is not too large. Consequently 
there is a correlation between the mass and condensate and we establish the dependence of 
the MS charm mass rh on the condensate (^G^) found at the point n = 10, = 0.98 x 4m^ 
on this diagonal. It is plotted in the fig The error of the experimental moments is about 
7%, arising mainly from the uncertainty in Fj/^^ee- But, since M„(Q^) ~ (4m^ + Q^)"", 
the mass error is of order 7/n%, i.e. is much smaller. For instance, at zero condensate 



m(m2) = 1.283 ±0.007 GeV for (^G'^) = (42) 



the error is purely experimental. The dependence plotted in fig ^ as well as the value (|4^) 
are weakly sensitive to particular choice of the QCD coupling «<, and the scale /i^. This is an 
obvious advantage of nonzero while the analysis at = leads to significantly higher 
error 



It is more difficult to find the restrictions on the mass and condensate separately. For this 
purpose one has to choose the point in (n, Q"^) plane which is 1) out of the diagonal, since no 
new information can be obtained from there, 2) not in the lower right corner (high n, low Q"^), 
where perturbative corrections as well as as corrections to the gluon condensate are large 
and 3) not in the upper left corner (low n, high Q^), where the continuum contribution to the 
experimental moments is uncontrollable. It turns out that if one considers the ratio of the 
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Figure 6: MS mass versus gluon condensate obtained from the ratio Mg/Mio above the 
diagonal. For more notations see Fig ^. 

moments (0)) ^^e mass-condensate dependence appears to be different in comparison to the 
fig|^. In particular, the results obtained from the ratio M10/M12 at = 4m^ and M15/M17 
at = Srh"^ are demonstrated in the left and right parts of the Fig ^ respectively. At the 
same figures the mass-condensate dependence, obtained from the moments Mio((5^ = 4m^) 
and Mi^lQ"^ = Srh"^) is also plotted. The error bands include both the experimental error of 
the ratio (|TTD and the uncertainty of (^Ol). Obviously the results, obtained outside the 
diagonal, are sensitive to the choice of as well as /i^. The small variation of /i^ slightly 
changes the acceptable region in the fig |, but if one takes /x^ few times lower, the region 
expands to the left significantly. 

The absolute limits of the MS charm quark mass and the gluon condensate can be 
determined from the fig |^: 

m(m2) = 1.275 ±0.015 GeV , {-^^^) = 0.009 ± 0.007 GeV^ (43) 

The restrictions on rh and the gluon condensate, obtained from other ratios of moments, 
agree with (^Bj), but are weaker (see Fig |^, where the ratio Mg/Mio is considered). The 
stability intervals in the moments, i.e. the intervals, where (|^) takes place within the errors, 
were found to be n = 8 — 13 at = 4m^ and n = 12 — 19 at = 8m^. 

As a check, the calculations were performed, where the a^-terms in the MS moments 
were omitted as well as the a^-corrections to the gluon condensate contribution (M„ = 
M^^'^) = 0). At = 4^2 ^g^g fQ^j^^ ^ ^ 1.266 GeV and {aG^) = 0.011 GeV^ from 
Mio and M^o/Mu while at = Srh"^ the result m = 1.263 GeV and (aG^) = 0.015 GeV^ 
was obtained from M15 and M15/M17. These values agree with ( ^3)) in the limit of errors. 
However, it is difficult to estimate the errors of the calculation, where terms are omitted, 
because of the uncertainty in the scale. 

Rc{s) in (13), or the expression for the moments are, in principle, independent on the 



normalization scale /i^. However, in fact, since we take into account only first 3 terms 
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in the tts-expansion in (|27|), such dependence takes place. Namely, when we change the 
normalization point in from fn^ to /i^ = + fn^ pi] ) with the help of eq (PU]), the values 
of the moments, defined by (pT]), are changed. As is clear from (p9|), the difference between 
the moments M„((5^) at the normalization points and + fn? increases with Q^. At 
used above, the difference is moderate and, if recalculated to (^^G'^), results in the error 
A (^G^) ~ 2 X 10~^ GeV^, much smaller, than the overall error in (^31) . However, going to 
the higher would be dangerous. In fact, while deriving (|39|), we expanded the running 
QCD coupling a(/x^) in a(m^). This expansion is valid if 

a/3oln^ < 1 (44) 

In particular for Q^/(4m^) = 3 the l.h.s. of this equation is ~ 0.5 and the neglected higher 
order terms could be significant. For this reason we avoid to use higher Q^, than it was done. 

Let us now turn the problem around and try to predict the width Tj/^^ee theoreticaly. 
In order to avoid the wrong circle argumentation we do not use the condensate value just 
obtained, but take the limitation (^G^) = 0.006 ± 0.012 GeV^ found in ^ from r-decay 
data. Then, the mass limits m = 1.28 — 1.33 GeV can be found from the moment ratios 
exhibited above, which do not depend on Tj/^^ee if the contributions of higher resonances 
is approximated by continuum (the accuracy of such approximation is about 3%). The 
substitution of these values of rh into the moments gives 

r*j);°:.,, = 4.9±0.8keV (45) 

in comparison with experimental value Fj/^^ge = 5.26 ±0.37keV. Such good coincidence of 
the theoretical prediction and experimental data is a very impressive demonstration of the 
QCD sum rules effectiveness. It must be stressed, that while obtaining (^) no additional 
input were used besides the condensate restriction taken from |T3 and the value of ^^(m^). 



7 D = 6 condensate influence 

The D = 4 gluon condensate {aG"^) is the leading term in the operator expansion series. 
The question arises, how the higher dimension condensate could change the results of our 
analysis. There is single D = 6 gluon condensate (g^G^). Its contribution to the polarization 
function (1^) can be parametrized as follows: 



The dimensionless function f^'^^\z) has been found in PD| : 

f'^'Kz) = ( - + - z + 4J2 - - J3 + - J, - - J, ] , (46) 

^ ' 12tx^z^ \lh 5 3 5 5 '' ^ ' 

where the integrals 

dx 

" " io [l-4^x(l-x)]'^ 
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Figure 7: MS mass versus gluon condensate obtained from the moments and ratios with 
account of (G^) condensate according to (|48|). 



can be calculated analytically. However the integral representation is convenient to express 
the result in terms of Gauss hypergeometric function, which can be easily differentiated in 
order to obtain the moments: 



(47) 



where 



1080 (4m2)"H 



r(n + i)r(n + 5) 

V{n + 1) T{n + 9/2) ' ' 



n + i, n + 5 
9/2 



n 



Am? 



and the constants C2 = 3, C3 = —7, C4 = —9. Significance of the condensate (g^G^) is 
determined by the ratio of the two terms in (^Tj). The numerical values of this ratio for 
different (n, Q^) are given in the last column of the Tables 1,2,3 in the Appendix. 

No reliable estimations of the (G^) condensate are available. There exists only the esti- 
mation, based on the dilute instanton gas model PT| : 



I „?> fahc/^a /~ib ric \ 
\9 J ^^lu^uX^Xfi/ 



4 127r^ /a,^2 

5 pI \ TT 



(4J 



where pc is effective instanton radius. The numerical value of pc is uncertain, even in the 
framework of the model: in ||3^ the value pc = 1/3 fm = 1.5GeV~^ was advocated, in 
the value pc = Ifm = 4.5GeV~^ was used. In the recent paper [^, based on the sensitive 
to gluon condensate sum rules, pc = 0.5 fm = 2.5 GeV~^ was suggested. 

The contribution of (g^G^) to MniQ"^) at a fixed n falls rapidly with growth of Q^. At 
Q2 = and n > 5 it comprises about 50% or more of the gluon condensate contribution at 
Pc = 0.5 fm. Even at (5^/(4m^) = 1 it is significant: the (negative) correction to the gluon 
condensate term is ~ 10% in Mio and ~ 30% in the ratio M10/M12. One gets more reliable 
results at Q^/(4m^) = 2. Here the corrections are: —7% for M15 and —18% for M15/M17. 
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These corrections leave the charm quark mass almost unchanged, but increase the gluon 
condensate and its error (compare Figs ^and|^). The account of {g^G^) contribution leads 
to the following restriction: 

(^^G^^ = 0.011 ± 0.009 GeV^ (49) 
Certainly, it relies upon the instanton gas model, that gives (^8]). 



8 About the attempts to sum up the Coulomb-hke cor- 
rections 

Sometimes when considering of the heavy quarkonia sum rules the Coulomb-like corrections 
are summed up ||3^-[0. The basic argumentation for such summation is that 

at = and high n only small quark velocities v < l/v^ essential and the prob- 
lem becomes nonrelativistic. So it is possible to perform the summation with the help of 
well known formulae of nonrelativistic quantum mechanics for |?/'(0)p in case of Coulomb 
interaction (see [0). 

This method was not used here for the following reasons: 

1. The basic idea of our approach is to calculate the moments of the polarization operator 
in QCD by applying the perturbation theory and OPE (l.h.s. of the sum rules) and to 
compare it with the r.h.s. of the sum rules, represented by the contribution of charmonium 
states (mainly by J/ip)- Therefore it is assumed, that the theoretical side of the sum rule 
is dual to experimental one, i.e. the same domains of coordinate and momentum spaces 
are of importance at both sides. But the charmonium states (particularly, J/ip) are by no 
means the Coulomb systems. A particular argument in favor of this statement is the ratio 
Fj/^^ee/rv'^ee ~ 2.4. If charmouia were nonrelativistic Coulomb system, T^^ee would be 
proportional to |?/'(0)p ~ ^/{nr + 1)^, and since is the first radial excitation with rir = 1, 



this ratio would be equal to 8 (see also ||38|). 

2. The heavy quark- ant iquark Coulomb interaction at large distances r > r^onf ~ 
1 GeV^^ is screened by gluon and light quark-antiquark clouds, resulting in string formation. 
Therefore the summation of Coulombic series makes sense only when the Coulomb radius 
rcoui is below rconf- (It must be taken in mind, that higher order terms in Coulombic series 
represent the contributions of large distances, r 3> rcoui-) For charmonia we have 

2 

rcoui - — ~ 4GeV 

mcCpas 

It is clear, that the necessary condition -Rcoui < -Rconf is badly violated for charmonia. This 
means that the summation of the Coulomb series in case of charmonium would be a wrong 
step. 

3. Our analysis is performed at Q^/Am"^ > 1. At large the Coulomb corrections are 
suppressed in comparison with = 0. It is easy to estimate the characteristic values of 
the quark velocities. At large n they are v ~ ^(l + QV4m2)/n. We are working along 
the diagonals of the Fig ^, well parametrized by the equation Q'^/Afh'^ ^ n/5 — 1. Here the 
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quark velocity v ~ ~ 0.45 is not small and not in the nonrelativistic domain, where 

the Coulomb corrections are large and legitimate. 

Nevertheless let us look on the expression of Rc, obtained after summation of the Coulomb 
corrections in the nonrelativistic theory It reads (to go from QED to QCD one has to 
replace a Cpas, Cp = 4/3): 

Rc Coul = ^ = -V [ I + - + + ... (50) 

ccoul 2 1-e-^ 2 V 2 12 720 / ^ ' 

where x = TrCpds/v. At f = 0.45 and as ~ 0.26 the first 3 terms in the expansion (|50|), 
accounted in our calculations, reproduce the exact value of -Rc, Coui with accuracy 1.6%. Such 
deviation leads to the error of the mass rh of order (1 — 2) x 10~^ GeV, which is completely 
negligible. In order to avoid misunderstanding, it must be mentioned, that the value of 
Rc Coul, computed by summing the Coulomb corrections in nonrelativistic theory has not 
too much in common with real physical situation. Numerically, at choosen values of the 
parameters, -RcCoui ~ 1-8, while the real value (both experimantal and in the perturbative 
QCD) is about 1.1. The goal of the argumens, presented above, was to demonstrate, that 
even in the case of Coulombic system our approach would have a good accuracy of calculation. 

At V = 0.45 the momentum transfer from quark to antiquark is Ap ~ 1 GeV. (This is 
typical domain for QCD sum rule validity.) In coordinate space it corresponds to Ar^g ~ 
1 GeV~ . Comparison with potential models demonstrates, that in this region the 



effective potential strongly differs from Coulombic one. 

4. Large compensation of various terms in the expression for the moments in MS scheme 
(see Fig |^) is not achieved, if only the Coulomb terms are taken into account. This means, 
that the terms of non-Coulombic origin are more important here, than Coulombic ones. 

For all these reasons we believe, that the summation of nonrelativistic Coulomb correc- 
tions is inadequate in the problem in view: it will not improve the accuracy of calculations, 
but would be misleading. 



9 Results and discussion 

The analysis of charmonium sum rules is performed within the framework of QCD at the 
next level of precission in comparison with famous treatment of this problem by Shifman, 
Vainstein and Zakharov [|I| . In the perturbation theory the terms of order were accounted 
as well as as corrections to the gluon condensate contribution, in OPE — the dimension 6 
operator G^. The method of the moments was exploited. The validity of the method was 
demonstrated for the MS mass of the charm quark, but not for the pole mass. The domain in 
the plane (n, Q'^) was found, where the three accounted terms in the perturbative series are 
well converging. It was shown, that the sum rules do not work at = 0, where the following 
4 requirements cannot be satisfied simultaneously: 1) convergence of the perturbation series, 
2) small as correction to the gluon condensate contribution, 3) small contribution of 
operator, 4) small contribution of higher resonances and continuum. Large allow also 
to suppress the Coulomb corrections. The most suitable values of for the sum rules are 
~ (1 — 2)47fi^ ~ 5 — 15 GeV^. The values of charmed quark MS and the gluon condensate 
were found by comparing the theoretical moments with experimental ones, saturated by 
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charmonium resonances (plus continuum). A strong correlation of the values m and (^G^) 
was established. This connection only weakly depends on Ug- Taking the value found in 



T9| from hadronic r-decay data 

as{ml) = 0.330 ±0.025, (51) 
the MS charm quark mass and the gluon condesate were determined 

m(m2) = 1.275 ±0.015 GeV , (—G^) = 0.009 ± 0.007 GeV^ (52) 



TT 



The error in ( ^2]) roughly comprises as 50% theoretical (uncertainly in as and the normal- 
ization scale) and 50% experimental (mainly the error of J/ip electronic decay width). The 
numbers in ( [5^ ) were obtained disregarding the contribution of operator. The account of 
G^ term, when {G'^) was taken using the dilute instanton gas model with pc = 0.5 fm, shifts 

m to 



—G^) = 0.01 1± 0.009 GeV^ (53) 

71 



The value ([53| ) may be compared with recently found |jT9[ limitation on the gluon condensate 
from hadronic r-decay data: 

^G^) = 0.006 ± 0.012 GeV^ (54) 

TT / 

Eqs ( |53D and (|5^ are compatible and obtained from independent sources. So, with some 
courage, we can average them and get 

= 0.0085 ± 0.0075 GeV^ (55) 

7T / av 

After such averaging we come back to (0). 

We can formulate our final conclusion about the gluon condensate value in such a way. 
The values of gluon condensate two times (or more) larger than the SVZ value (H) are 
certainly excluded. Unfortunately our analysis does not allow to exclude zero values of 
the gluon condensate. In this aspect the improvement of the experimental precission of 
J/ijj —>■ e^e~ width would be helpfuU. Based on the condensate limitation (Q) and the 
value of as (^iD, the J/iIj electronic decay width Fj/^^ee was predicted theoretically: 

r*j);°:.,, = 4.9±0.8keV (56) 

in comparison with the experimental value 5.26 ± 0.37 keV. Such a good coincidense ones 
more demonstrates the effectiveness of QCD sum rule approach. 
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Appendix: Numerical values of the moments 

We list here the numerical values of the perturbative moments M*^°'^'^\ condensate contribu- 
tion M(<^'0'1) in MS scheme computed by (|37|) and (G^) condensate contribution M^"^^^ (|47|) 
for n = 1..20 and Q'^/ (4m^) = 0, 1, 2. For dimensionfull values we put 4rn^ = 1 here, so that 
the leading term M(°) and the ratios M^^'^y M^^\ M^^^y M'^^^ should be divided by {Am^Y 
and (4m^)^, (4m^)^ respectively for a particular mass m. 
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2 
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32.21 


6 


7.9565 
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10- 


-2 


-1.092 


0.213 


-780.3 
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-2.375 
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-3011. 
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16 
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-2 
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86.51 


-1.103 X 10^ 


-24.81 
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17 
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-2 


-18.49 
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-1.308 X 10^ 
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18 
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-1.537 X 10^ 
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X 
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-2 


-22.13 
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-1.791 X 10^ 


-30.98 


1.268 X 10^ 
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-2 


-23.98 


205.2 


-2.072 X 10^ 


-33.07 


1.614 X 10^ 



Table 1: Moments at = 
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10- 


-3 


2.34 
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12 


9.5483 X 


10- 


-6 


-1.167 


-1.656 


-700.3 


-3.346 


38.5 


13 


4.2743 X 


10- 


-6 


-1.815 


-2.732 


-867.3 


-4.205 


56.96 


14 


1.9283 X 


10- 


-6 


-2.486 


-3.508 


-1059. 


-5.078 


81.75 


15 


8.7574 X 


10- 


'7 


-3.178 


-3.93 


-1276. 


-5.962 


114.4 


16 


4.0007 X 


10- 


-7 


-3.89 


-3.948 


-1521. 


-6.858 


156.4 


17 


1.8372 X 


10- 


-7 


-4.62 


-3.518 


-1795. 


-7.764 


209.9 


18 


8.4756 X 


10- 


-8 


-5.365 


-2.599 


-2101. 


-8.68 


277. 


19 


3.9264 X 


10- 


-8 


-6.126 


-1.154 


-2439. 


-9.605 


360. 


20 


1.8257 X 


10- 


-8 


-6.9 


0.85 


-2811. 


-10.54 


461.7 



Table 2: Moments at = Afh'^ 
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n 














1 


0.3005 


2.073 


-2.016 


-1.098 


5.002 


-8.936 X 10-3 


2 


4.6172 X 10-2 


2.531 


2.453 


-3.825 


4.762 


-0.03 


3 


9.589 X 10-=^ 


2.734 


5.208 


-8.691 


4.468 


-0.064 


4 


2.2613 X 10-^ 


2.792 


6.909 


-16.2 


4.131 


-0.105 


5 


5.7274 X 10"^ 


2.753 


7.853 


-26.84 


3.761 


-0.14 


6 


1.5191 X 10"^ 


2.643 


8.229 


-41.11 


3.364 


-0.147 


7 


4.1621 X 10"^ 


2.478 


8.167 


-59.5 


2.942 


-0.094 


8 


1.1682 X 10-^ 


2.27 


7.769 


-82.51 


2.501 


0.066 


9 


3.3404 X 10-6 


2.025 


7.113 


-110.6 


2.042 


0.393 


10 


9.6957 X 10"^ 


1.749 


6.264 


-144.3 


1.568 


0.963 


11 


2.8488 X 10"^ 


1.447 


5.276 


-184.1 


1.08 


1.871 


12 


8.4559 X 10-s 


1.122 


4.195 


-230.4 


0.579 


3.233 


13 


2.5317 X 10-^ 


0.776 


3.061 


-283.9 


0.068 


5.186 


14 


7.6361 X 10-9 


0.412 


1.909 


-344.8 


-0.455 


7.89 


15 


2.3181 X 10-9 


0.031 


0.77 


-413.8 


-0.986 


11.53 


16 


7.0768 X 10-^° 


-0.364 


-0.33 


-491.3 


-1.527 


16.33 


17 


2.1713 X 10-1° 


-0.773 


-1.365 


-577.8 


-2.075 


22.52 


18 


6.6914 X 10-11 


-1.195 


-2.313 


-673.9 


-2.631 


30.38 


19 


2.0704 X 10-11 


-1.628 


-3.153 


-779.9 


-3.194 


40.21 


20 


6.429 X 10-12 


-2.072 


-3.867 


-896.4 


-3.764 


52.36 



Table 3: Moments at ^ 3^2 



24 



